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Calculating e'/e using HYP staggered fermions* 
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We present preliminary results for e'/e calculated using HYP staggered fermions in the quenched approximation. 
We compare different choices of quenched penguin operators. 



There are a number of advantages to using 
staggered fermions for calculating weak matrix 
elements relevant to CP violation in the neutral 
kaon system. They retain sufficient chiral sym- 
metry to protect the weak operators from mix- 
ing with lower dimensional operators. Simulat- 
ing them is considerably cheaper than domain 
wall and overlap fermions. On the other hand 
unimproved staggered fermions have a number of 
drawbacks: 1-loop corrections to operator renor- 
malization are large (« 100%), scaling violations, 
even though they begin at order a 2 , are large as is 
the breaking of the SU(4) taste symmetry. Some 
of these limitations can be alleviated by improv- 
ing staggered fermions using fat links. 

Explicit 1-loop calculations for bilinear opera- 
tors show that taste symmetry breaking and the 
renormalization corrections can be reduced by an 
order of magnitude by improving the lattice ac- 
tion. The greatest reduction is observed for Fat7 
and mean field improved HYP/Fat7 pQ. Sim- 
ilarly, the complete set of 1-loop matching co- 
efficients for the four-fermion operators relevant 
to CP violation also have small corrections at 
1/a m 2 GeV |1I5| . Lastly, numerical simulations 
show that HYP smeared links reduce taste sym- 
metry breaking in the pion multiplet 2 . In light 
of these results, we have chosen to use HYP/Fat7 



* Presented by W. Lee. Research supported in part 
by BK21, by Interdisciplinary Research Grant of Seoul 
National University, by KOSEF contract R01-2003-000- 
10229-0, by US-DOE contracts DE-FG03-96ER40956, W- 
7405-ENG-36, KA-04-01010-E161, and NERSC. 



staggered fermions to calculate matrix elements 
relevant to e'/e. 

In the standard model, e'/e is: 



e'/e = 

p(l/2) = 
p(3/2) = 



im(v; s v td ) 
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where Vij are elements of the CKM matrix, 
to = 1/22.2 quantifies the AI — 1/2 rule and 
fl v +r]' = 0.060(77) represents the isospin break- 
ing effect 0. P^ 1 / 2 ) (P( 3 / 2 )) are the AI = 1/2 
(AI — 3/2) contribution to e'/e and the yi(n) are 
Wilson coefficients in the operator product ex- 
pansion, given in JJJ. The sum over operators Oi 
includes the QCD (i = 3,4,5,6) and electroweak 
(i = 7,8,9,10) penguin operators. The current- 
current operators (i = 1,2) do not contribute to 
CP violation. On the lattice we calculate the ma- 
trix elements (ir\Oi(n)\K) and {0\Oi(/j,)\K ) and 
relate them to (Oi)i(^) = (Trnj\Oi(^i)\K) using 
chiral perturbation theory at the leading order. 
Here the subscript / refers to the isospin of the 
final two pion state. 

Matrix elements are calculated on 218 
quenched gauge configurations at [3 — 6.0 



64 lattices. Results arc converted 
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Figure 1. Contributions to p( x / 2 ), using linear 
chiral fits. Empty (shaded) columns correspond 
to the STD (GP) methods for (Q 5 , 6 ). 

into the continuum NDR scheme (matching at 
q* = 1/a = 1.95 GeV), evolved to the charm 
quark mass m c = 1.3 GeV (using running with 
Nf — 3) [7| and then combined with the Wilson 
coefficients yi(m c ) to obtain P^ 1 / 2 ) and p( 3 / 2 ). 

We compare two ways of transcribing opera- 
tors O5 and Oe to the quenched theory. In the 
first (the standard [STD] method), the operators 
belong to the singlet representation of SU(3)r, 
while in the second (the Golterman-Pallante [GP] 
method |H]) they belong to the singlet represen- 
tation of the graded group SU (3\3)r. They dif- 
fer because a singlet under SU(3)r is a combina- 
tion of singlet and non-singlet representations of 
the SU (3\3)r. Since the dominant contribution 

tO P(V2) 

comes from Oq, the difference between 
these two methods can and, as we show, does have 
a significant impact on e'/e. 

A second important issue in the analysis is the 
function used to extrapolate the data to the chiral 
limit. First, we take an appropriate ratio to ob- 
tain quantities which does not vanish in the chiral 
limit. We would then like to fit the data to the 
general form 

ci + c 2 m 2 K + c 3 m 2 K \og(m 2 K ) + CA{m 2 K ) 2 , (5) 

except that we lack measurements over a suffi- 
cient range and number of masses to simultane- 
ously determine all the c; parameters. Thus, we 
either exclude the logarithm (03 = 0: "quadratic 
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Figure 2. Contributions to p' 1 / 2 ), using linear 
(left columns) and quadratic (right columns) chi- 
ral fits, and the GP method for (Qsfi). 

fit"), or the quartic term (04 = 0: "log fit"), or 
both (03 = C4 = 0: "linear fit"). In general, 
our data cannot distinguish between the fits. In 
most cases there is no useful theoretical guidance 
for the parameter of the fits, and we compare all 
three. 

Results for P^ 1 / 2 ) using the both STD and GP 
transcriptions of (O5} and (Og) and linear extrap- 
olation to the chiral limit are shown in Fig. ^ 
The last column in Fig. ^ is the sum of individ- 
ual contributions. The dominant contributions 
are from {Oq) and (O4), which partially cancel. 
Note that the evolution from q* to m c mixes O5.6 
with other operators, so all contributions depend 
on the choice of the operator. Clearly the choice 
of operator makes a large difference to the final 
result. 

Figure [3 illustrates the dependence of the re- 
sults on the fitting function used in the chiral ex- 
trapolation. The errors increase on going from 
linear to quadratic to log fit so much so that the 
difference in the result is only about one stan- 
dard deviation. The dominant source of the un- 
certainty is the value and error on the point at 
the lightest quark mass. This is illustrated in 
Fig. [3] which compares linear and quadratic chiral 
extrapolations for the ratio {Qe) / {m 2 K f 2 ), which 
should be non-vanishing in the chiral limit. The 
log fit, which is not shown to maintain clarity, ex- 
trapolates to —30(12). So, on the one hand, one 




Figure 3. Linear and quadratic fits to 
(ir\Q 6 \K)/{m 2 K f 2 ) (GP method). 

needs masses at least this light (w m s /5) to do 
reasonable chiral extrapolations, but on the other 
hand the statistical errors and possibly finite vol- 
ume effects grow rapidly with 1/m. 

Lastly, in Fig. 21 we show the contributions to 
p(3/2) usm g linear fits. Here there are no chiral 
logarithms in the quenched theory [S] and there 
is no ambiguity in the operator definition. The 
contribution from (Og) dominates, as expected 
from a large N c analysis. 

Using Eq. Q our preliminary estimates for e'/e 
for the STD operators are 



(STD) 



-2.4(34) x 1CT 4 (linear fit) 
-12.8(114) x 1(T 4 (quad, fit) 
-20.1(175) x 10~ 4 (log. fit) 



These results are consistent with previous esti- 
mates obtained using domain wall fermions (and 
STD operators) PTU| . On the other hand the GP 
method gives a more positive estimate for (Oe), 
and consequently e'/e: 



(GP) 



+3.2(29) x 1CT 4 
+8.8(93) x 10~ 4 
+13.1(142) x 10~ 4 



(linear fit) 
(quad, fit) 
(log. fit) 



The very large errors in these results are caused 
by the cancellation between the P^ 1 / 2 ) and p( 3 / 2 ) 
contributions. Nevertheless, it appears that the 
difference between GP and STD, which is a 2 
sigma effect, is significant. To clarify this point, 
and to improve the chiral extrapolations, we are 
extending the current run to 400 configurations 



Figure 4. Contributions to p( 3 / 2 ) . 

with 4 additional quark masses. Ultimately, of 
course, one needs to do partially quenched simula- 
tions where the operator ambiguity is less impor- 
tant and the coefficients of the chiral logarithms 
are known. 

This calculation is being done on the Columbia 
QCDSP supercomputer. We thank N. Christ, 
C. Jung, C. Kim, G. Liu, R. Mawhinney and 
L. Wu for their support on the staggered e'/e 
project. 
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